Brownian motion with respect to time-changing 
Riemannian metrics, applications to Ricci flow 



K.A. Coulibaly 



Abstract 

We generalize Brownian motion on a Riemannian manifold to the case of 
a family of metrics which depends on time. Such questions are natural for 
equations like the heat equation with respect to time dependent Laplacians 
(inhomogeneous diffusions). In this paper we are in particular interested in 
the Ricci flow which provides an intrinsic family of time dependent metrics. 
We give a notion of parallel transport along this Brownian motion, and es- 
tablish a generalization of the Dohrn-Guerra or damped parallel transport, 
Bismut integration by part formulas, and gradient estimate formulas. One of 
our main results is a characterization of the Ricci flow in terms of the damped 
parallel transport. At the end of the paper we give a canonical definition of 
the damped parallel transport in terms of stochastic flows, and derive an in- 
trinsic martingale which may provide information about singularities of the 
flow. 



g(£)-Brownian motion 



Let M be a compact connected n-dimensional manifold which carries a family of 
time-dependent Riemannian metrics g(t). In this section we will give a generaliza- 
tion of the well known Brownian motion on M which will depend on the family of 
metrics. In other words, it will depend on the deformation of the manifold. Such 
family of metrics will naturally come from geometric flows like mean curvature flow 
or Ricci flow. The compactness assumption for the manifold is not essential. Let 
V* be the Levi-Civita connection associated to the metric g(t), A t the associated 
Laplace-Beltrami operator. Let also (fi, (J-t)t>o, J~, P) be a complete probability 
space endowed with a filtration (-7~i)t>o satisfying ordinary assumptions like right 
continuity and W be a IR n -valued Brownian motion for this probability space. 

Definition 1.1 Let us take (Q, (J r t)t>o, J 7 , P) and a C 1 ' 2 -family g(t)te[o,T[ of met- 
rics over M. An M -valued process X(x) defined on f2 x [0,T[ is called a g(t)- 
Brownian motion in M started at x G M if X(x) is continuous, adapted, and if 
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for every smooth function f , 

f(X s (x)) - f(x) - ± jT AJ(X t (x)) dt 

is a local martingale. 

We shall prove existence of this inhomogeneous diffusion and give a notion of 
parallel transport along this process. 

Let (e^jgji..^ be an orthonormal basis of M™, JF(M) the frame bundle over M, ti 
the projection to M. For any u G JF(M), let Li(t, u) = /i*(uej) be the V* horizontal 
lift of uti and Li(t) the associated vector field. Further let V a> p be the canonical 
basis of vertical vector fields over T{M) defined by V a ^{u) = Dl u (E a ^) where E a> p 
is the canonical basis of A4 n (M.) and where 

l u : GL n (R) -> JT(M) 

is the left multiplication. Finally let (O(M), g(t)) be the g{t) orthonormal frame 
bundle. 

Proposition 1.2 Assume that g(t)t<=[o,T{ is a C 1,2 (t,x) -family of metrics over M, 
and 

A : [0,T[ x T(M) -> M n (R) 

(t,U) i-> (A ajP (t,U)) ajl3 

is locally Lipschitz in U unformly in all compact oft. Consider the Stratonovich 
differential equation in JF(M); 

f *dU t = Y:tiHt,Ut)*dW i + j: a ^A a>l3 (t,U t )V a>l3 (U t )dt 

\U e F(M) such that U E (O(M),g(0)). [ ' ' 

Then there is a unique symmetric choice for A such that U t G (O(M), g(t)). More- 
over: 

A{t i U) = ~d 1 G(t,U) i 
where {diG(t,U)) it j = (Ue u Uej) dtg{t) . 

Proof: Let us begin with curves. Let / be a real interval, 7r : TM — > M the 
projection, V and C in C 1 (/,TM), two curves such that 

x(t) := 7r(V(t)) = n(C(t)), for all t € / 

We want to compute: 
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We write dig(t,x) for d s g(s,x) evaluated at t. Let us express the metric g(t) in 
a coordinate system; without loss of generality we can differentiate at time 0. Let 
(x 1 , ...,x n ) be a coordinate system at the point x(0), in which we have: 

V(t) = v\t)d x , 

c(t) = c\t)d x i 

g(t, x(t)) = gi,j(t, x(t))dx l (g) dx j 
In these local coordinates we get: 

4 (nt),C(t)) g{tAt)) = | g h ^x{t))v\t)c>{t) 
ai | t=0 at | (=0 

= (9 ia ,(0,x)^(0)tf(0) + {g t ^x{t))v\ty{t)) 

at U_n 



^■(0,^(0)^(0) + ( V° (o) V(0), C(0) 

I g(0,x(0)) 

In order to compute the g{t) norm of a tangent valued process we will use what 
Malliavin calls "the transfer principle", as explained in [T3].[T2]. 

Recall the equivalence between a given connection on a manifold M and a 
splitting on TTM, i.e. TTM = H V TTM © VTTM p]. We have a bijection: 

V, : T n(v) M — > V V TTM 

u ^ s(« + *«)lt=o- 

For X, y G T(TM) we have: 

v JC y(x) = v^„ ) ((dy(x)(x(x))) w ), 

where (.) v is the projection of a vector in TTM onto the vertical subspace VTTM 
parallely to H^TTM. 

For a T(M)-valued process T t , we define: 

D s '%=(V Tt )- 1 {{*dT t ) v ' t ), (1.2) 

where (.)"'' is defined as before but for the connection V*. The above generalization 
makes sense for a tangent valued process coming from a Stratonovich equation like 
U t Si, where U t is a solution of the Stratonovich differential equation (II. lft . 
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For the solution U t of (jl.lD we get 



d (( U t e u u te 3 ) g{tMUt)) ) = (UteuUte^g^^dt (1.3) 

+ <D**l7 tei) ^•> g( , 7r(f7t)) + (U tei , D s '%e,) g{tM ^f) 

We would like to find a symmetric A such that the left hand side of the above 
equation vanishes for all time (i.e. U t G (O(M), g(t))). Denote by ev £i : JF(M) — > 
TM the ordinary evaluation, and dev £i : TT{M) —> TTM its differential. 
It is easy to see that dev ei sends VTT[M) to VTTM and sends H vh TT[M) to 
H^TTM. We obtain: 

r( 

D s ' t U t e i = ^A a ^U t )U t e a dt. (1.5) 

a=l 

For simplicity, we take for notation: (diG(t, U)) id = (Uei, Uej) dtg ^ and 

(G(t, E0)y = (f/e i; f/ ej ) fl(t) . 

It is now easy to find the condition for A: 

(G(t, U t )A(t, U t )) hl + (G(t, U t )A(t, U t )) id = -(d 1 G(t, U t )) id (1.6) 

Given orthogonality G(t,U t ) = Id and so by (jl.6l) A differs from — by skew 
symmetric matrice, therefore will be equal to it if we demand symmetry. Conversely 
if A = —\d\G then by ( 11 .31) and equation (11.21) we see G(t, U t ) = Id. 
□ 

Remark : The SDE in proposition 11.21 does not explode because on any 
compact time interval all coefficients and their derivatives up to order 2 in space 
and order 1 in time are bounded. 

Remark : The condition of symmetry is linked to a good definition of parallel 
transport with moving metrics in some sense. 

To see where the condition of symmetry comes from we may observe what 
happens in the constant metric case. It is easy to see that the usual definition 
of parallel transport along a semi-martingale which depends on the vanishing of 
the Stratonovich integral of connection form, is equivalent to isometry and the 
symmetry condition for the drift in the following SDE in !F{M): 

dUt = Eti HUt) * dW l + A(U t ) atP V a AUt) dt 
U e(O(M),g) 

U t E (0(M),g) (isometry) 

-) a ,f3 £ S{n) (vertical evolution). 
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Remark : Isometry of U t forces A to be skeew symmetric by (11.61) . the symmetry 
of A give A = 0. We get the usual stochastic differential equation of the parallel 
transport in constante metric case. 

The next proposition is a direct adaptation of a proposition in [T5], page 42; 
hence the proof is omitted. 

Proposition 1.3 Let a £ Y(T*M) and F a : F{M) -> R d , F^(u) = a w{u) (uei) its 
scalarization. Then, for all A £ T(TM), 

(V j4 tt) 7r ( M )(ue i ) = h(A w{u) )F l a . 

Consequently, for all u £ JF(M), 

(Vf UfUu^ue,) = h 9 ®(A <u) )F l df 

and for / £ C°°(M), 

Li(t)(foTr){u) = d{foTr)Li{t,u) 
= Fjf(u). 

Hence we have the formula: 



(v£?#)( 



tte. 



= W^dfiue^uej). 
Proposition 1.4 Take x £ M and the SDE %n T{M): 

\ *du t = J2ti L i(ti u t) * dw* - ^dMt, u t ) a , p v a> p{u t ) dt 

\ U £ F{M) such that U £ (O x (M), g(0)). [ > 

Then X t (x) = 7r(U t ) is a g(t)-Brownian motion, which we note g(t)-BM(x). 

Proof : For / £ C°°(M), 
d(fo7ro U t ) = J27=i L i(t)(f ° n)(U t ) * dW l 

= EILi HW o nWJdw* + i E- i= i hWAW ° ^)dw^ 

= lEti^dfiUte^U^dt 
d = lA t f(7roU t )dt. 

The last equality comes from the fact that U t £ (O(M), g(t)). □ 
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Remark : Recall that in the compact case the lifetime of equation fll.7j) is 
deterministic and the same as the lifetime of the metrics family. 



Let Ut be the solution of (11.71) . We will write //o,t = U t o Uq 1 the g(t) parallel 
transport over a g(t)-Brownian motion ( we call it parallel transport because it is 
a natural extention of the usual parallel transport in the constante metric case). 
As usual it is an isometry: 

// 0!t :(T Xo M,g(0))^(T Xt M,g(t)). 

We also get a development formula. Take an orthonormal basis (vi,...,v n ) of 
(Tx M, g(0)), and X t (x) a g(t)-Brownian motion of proposition [L~4l then 

*dX t (x) = // ,tVi * dW\. 

For / G C 2 (M) we get the Ito formula: 

df(X t (x)) = (W, /MtdW* + ~A t (f)(X t (x)) dt. (1.8) 

We will now give examples of g(t)-Brownian motion. Let (S n ,g(0)) be a sphere 
and the solution of the Ricci flow: -§^g(t) = — 2 Ricj that is g(t) = (l — 2(n — l)t)g(0) 
with explosion time T c = 2 ( n -x) ■ We will use the fact that all metrics are conformal 
to the initial metric to express the g(t)-Brownian motion in terms of the g(0)- 
Brownian motion. Let / G C 2 (S n ), X t (x) be a g(t)-Brownian motion starting 
at x G S n , B t some real-valued Brownian motion, and M t (x) a S n valued g(0)- 
Brownian motion. Then: 

df(X t (x)) =|| V*/(**(*)) IU) dBt + \ { j^—^t ) A of(X t (x))dt. 
We have: 

1 



Let 



then 



V</ l-2(n-l)t 11 "° • 



T(t)= / 1- 2(1-1)8^ 



ln(l-2(n-l)t) > T - 1(t) _e^ 



-2(n-l) -2(n-l) 
We have the equality in law: 

(X.(x)) = (B r( .)(x)). 
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We have a similar result for the hyperbolic case: Let (H n (— 1), g(0)) be the 
hyperbolic space with constant curvature —1. Then g{t) = (1 + 2{n — l)t)g(O) is 
the solution of the Ricci flow. Let X t (x) be a g(t)-Brownian motion starting at 
x G S^and M t (x) an H n -va\ued g(0)-Brownian motion. Then: 

/** 1 

r(t) = / -, r— ds, 

1 ; J l + 2(n-l)s 

and in law: 

(X.(x))±(M tQ (x)). 

Let us look at what happens for some limit of the Ricci flow, the so called 
Hamilton cigar manifold ([5]). Let on M 2 , g(0,x) = j^plp 9c&n be the Hamilton 
cigar , where || . || is the Euclidean norm. Then the solution to the Ricci flow 
is given by g(t,x) = -£rp^g{0, x). Let / G C 2 (IR 2 ), X t (x) be a g(t)-Brownian 
motion starting at x G IR 2 , B t a real- valued Brownian motion, and M t (x) some M 2 
valued g(0)-Brownian motion. Then: 

1 P 4t_|_ II vJr) II 2 

df(X t (x)) =|| VV(X t (x)) |U) dft + - i 2 A of(X t (x)) dt. 



We have: 



e 4 *+ II x 11 



WW = ; ; , v°/(x) 



1+ X 



We set: 



Then in law: 



p4£ _i_ II ~, 1 1 2 

v7(*) 11?= II v°/(x) ||g, 

p 4f , || 112 

A ^ = ITT i 
1+ x M 



(X. (x)) = (B r( .)(x)) 



Remark: If is a g(t)-Brownian motion associated to a Ricci flow 

started at g(Q) then X f / C (x) is a cg(t/c)-Brownian motion associated to a Ricci 
flow started at cg(0) so it is compatible with the blow up. 
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2 Local expression, evolution equation for the den- 
sity, conjugate heat equation 

We begin this section by expressing a g(t)-Brownian motion in local coordinates. 

Proposition 2.1 Let x G M, (x 1 , ...,x n ) be local coordinates around x, and X t (x) 
a g(t)-Brownian motion. Before the exit time of the domain of coordinates, we 
have: 

dXi(x) = yJg-^t,X t (x)) itj dB^ - l -g k > l r kl {t,X t {x))dt 

where we denote ^g _1 (t, X t (x))ij the unique positive square root of the inverse to 
the matrix (g(t, d x i, d x j))ij; here T kl (t, X t (x)) are the Christojfel symbols associated 
to V 9(t \ and B i are n independent Brownian motion. 

Proof: From the Ito equation 11.81 we get: 

dXl(x) = (VV, //o,tVi) gi t)dW l + ^ t x\X t {x))dt, 

where (vx, v n ) is a g(0)-orthogonal basis of T X M. By the usual expression of the 
Laplacian in coordinates: 

A t x*(X t (x)) = -g l > k r kl (t,X t (x)), 

and the gradient expression of the coordinates functions: 



we have: 



d , 1 



dXl(x) = g(t)^(^JM m dW l --g^Ti l (t,X t (x))M 

I 

dxj 



hit) 





£ VsTiVsPo-, /ktvi) m dw l - l -g l > k r u {t, x t (x)) dt 



Jg-(^dB m - ^g l > k r kl (t,X t (x))dt. 



where dB m = (^/^(t)" 1 '-?^-, //o,t v i)g(t)dW l . By the isometry property of the paral- 
lel transport and Levy's theorem B = (B 1 , B n ) is a Brownian motion in R n . 
□ 
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Remark : The above equation is similar to the equation in the fixed metric 

case. 



Now we shall study the evolution equation for the density of the law of the 
g(t)-Brownian motion. Let X t {x) be a p(i)-BM(x), and d\i t the Lebesgue measure 
over (M,g(t)). X t (x) is a diffusion with generator A t , we have the smoothness of 
the density (e.g. [22]). Let h x (t,y) £ C°°(]0,T[xM) such that: 

X t (x) = h x (t,y)diit(y),t>0 
X (x) = S x . 

By the continuity of X t {x) and the dominated convergence theorem we get the 
convergence in law: 

lim X t (x) = 5 X . 

We write in a local chart the expression of dfi t in terms of dfio, i.e., 



dfH = ^^'^K /detjg^midx 1 A dx 2 A ... A dx n \ 



VdetQ7 M (0)) 
and we note: 

Ht(dy) = ip(t,y)fi (dy). 

Proposition 2.2 

j t (h%t, y)) + h%t, y) Tr QgT^, 2/))|<?(*, 2/)) = ^A g(t) ^(t, y) 
c 

lim h x (t,y)dfi t = S x . 
Proof: For / £ C°°(M), t > 0, by definition of X t (x) we have: 



E[f(X t (x))] - f(x) = m 



J*A g{s) f(X s (x))ds 



f t E[f(X t (x))} = m[A g{t) f(X t (x))], 



dt 

i.e.: 

if M hx ( t >y)f(y)vt( d y) = U M A 9(t)f(y)h x (t,y)Mdy) 

The last equality comes from Green's theorem and the compactness of the manifold. 
By changing fJ>t(dy) = ip(t,y)fio(dy) in the left hand side, we have: 

f(v)l(h'(t,yMt,y))to(dy) = \ [ f(y)(A g{t) h x (t,ymt,y))^(dy) 
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so: 

j t (h*(t, y)i>(t, y)) = \(A g{t) h*(t, y)m, y)) (2.1) 
We also have by determinant differentiation: 

= ^(t,y)tr^g" 1 (t,y)j t g(t,y)y 

The part Tr g' 1 ^, y)ja9{t, vi^j is intrinsic, it does not depend on the choice 

of the chart. Hence ( 12.11) gives the following inhomogeneous reaction-diffusion 
equation: 

f t (h*(t,y)) + h*(t,y)tr(^g-\t,y)j t g(t,y^ = ±A g{t) h*(t,y). 

□ 

We will give as example the evolution equation of the density in the case where 
the family of metrics comes from the forward (and resp. backward) Ricci flow. 
From now Ricci flow will mean (probabilistic convention): 

|s M = -Ric M . (2.2) 

(respectively) 

Tt9 l j = R< j - (2-3) 
Remark: Hamilton in [H], and later DeTurck in [7] have shown existence 
in small times of such flow. In this section we don't care about the real existence 
time. 

For x G M, we will denote by S(t, x) the scalar curvature at the point x for the 
metric g{t). 

Corollary 2.3 For the backward Ricci flow (12.31) . we have: 

f t (h x (t,y)) + h*(t,y)S(t,y) = U g{t) h*(t,y) 
c 

lim h x (t,y)d/j t = 5 X . 
For the forward Ricci flow (12.21) . we have: 

f t (h x (t,y)) - \h*{t,y)S{t,y) = ±A g{t) h*(t,y) 
c 

lim h x (t,y)dnt = 5 X . 
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Remark : These equations are conservative. This is not the case for the ordi- 
nary heat equation with time depending Laplacian i.e. A 9 ( t ). They are conjugate 
heat equations which are well known in the Ricci flow theory (e.g. [24]). 



3 Damped parallel transport, and Bismut formula 
for Ricci flow, applications to Ricci flow for sur- 
faces 

In this section, we will be interested in the heat equation under the Ricci flow. The 
principal fact is that under forward Ricci flow, the damped parallel transport or 
Dohrn-Guerra transport is the parallel transport defined before. The deformation 
of geometry under the Ricci flow compensates the deformation of the parallel trans- 
port (i.e. the Ricci term in the usual formula for the damped parallel transport 
in constant metric case see ([9], [23], [El)). The isometry property of the damped 
parallel transport turns out to be an advantage for computations. In particular, for 
gradient estimate formulas, everything looks like in the case of a Ricci flat manifold 
with constant metric. We begin with a general result independent of the fact that 
the flow is a Ricci flow. Let g{t)[o t T c [ be a C 1,2 family of metrics, and consider the 
heat equation: 

/ d t f(t,x) = ±A t f{t,x) f . 

\ /(0,x) = / (x), 

where fo is a function over M. We suppose that the solution of (13.11 ) exists until 
T c . For T < T c , let Xf be a g(T — t)-Brownian motion, //o t the associated parallel 
transport. 

Definition 3.1 We define the damped parallel transport W^ t as the solution of: 
*d{{//l t )-\Wl t )) = -^(//y^Ric^) -d t (g(T - t)))* 9 V-HKt) dt 

with 

Wl t : T X M — > T X T {X) M } < = ld TxM . 
Theorem 3.2 For every solution f(t, .) of (13.11) . and for all v e T X M, 

df(T-t,.) X T (x) (Wl t v) 

is a local martingale. 
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Proof: Recall the equation of a parallel transport over the g(T — t)-Brownian 
motion Xj{x): 

*dUj = Eti HT - t, Uj) * dW* - \d t {g{T - t))(U?e a , U? e p )V a , p (U?) dt 
Uq e (O x (M),g(T)). 

(3.2) 

For / G C°°(M), its scalarization: 

df : F{M) — ► R n 

U i— ► (df(U ei ),...,df(Ue n )), 

yields the following formula in R n : 

df(T-t,.) xnx) (wl t v) = {df{T -t^^ujy^ltv)^ 

for every v e T X M. To recall the notation let: 

ev ei : F(M) — > TM 
[/ i — > Ua 

and recall that Uj , solution of (13.21) . is a diffusion associated to the generator 
\^_ t - l -8 t {g{T ~ *))(ev ei (.), cv, ,(.))l ;.,(.) 

where A^_ t is the horizontal Laplacian in A4, associated to the metric g(T — t). 
In the Ito sense, we get: 

d(df(T - 1, .) xTix) (wl t )v) = d(df(T - 1, uj), (ufy'wlv)^ 

^ (-(j/f)( T - *. -Wt)dt + [\^Jf{T - t, .) 

- \d t {g{T ~ t))(ev« , ev ej .)V iy j(.)df(T - t, .)}{Uj) dt, (U^W^v)^ 



+ <(#(T - t, t/f )), (^^((filWa),)^ 

f ~(j t df)(T - t, .)((W*» dt + (\\^Jf{T - t, .) 

- \d t {g{T ~ t))(ev« , ev e ,)V^(.)df(T - t, .)](f/ t T ) dt, (f/fl^Wj.^r 

- ^<(*(T - f, Uf)l {U^Y\//l t )-\Ri Cg{T - t) -d t (g(T - t)))*^(Wl t )vdt) Rn . 

We shall make separate computations for each term in the previous equality. Using 
the well known formula (e.g. [15], page 193) 

A H df = Adf, 
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we first note that: 

{\^Jf{T - t, .){UJI {UjY^l t vdt) v 

= |(A T _ t d/(T - t, .)(C/f), (U^Wl^n dt 
= 1 1 A T „ t df(T-t,.)(Wl t v) dt, 

By definition: 

Vijdf{u) = ^odfHld+tEij)) 

= ||*=o(t^(u(Id+t^) e -))»=i..» 
= (df(u5 s i e j )) s=1 .. n 

= (0, 0, df(uej), 0, 0) z-th position, 

so that: 

E« - <))(ev« ■, ev ei - *, -){UJ) dt 

= Eij dt(g(T - t)){Ufe u Vje 3 )df{Uje,)e t dt 

= ((V T -*/(T - t, .), Ei ®t(g(T - t))(U?ei, UfeM^T-t dt) i=l .. n 

= (df(T - t, d t (g(T - f))# r -*(0?e«)) d*) i=1 .. n . 

Then 

d(df(T-t,.) xT(x) ((Wl t )v)) 
d i?-idf(T-t,.)((Wl t v)dt 

-U(df(T - t, d t (g(T - t))# T -\Ure^)) t=1 .. n , (U^W^v)^ dt 
+|A r - t d/(T-t,.)(Wj t i;)dt 

-|((d?(T - t, C/f )), (^^(//Str^^cr-*) - ft(y(T - *)) # ^(Wj> dt> R . 
By the fact that Uj is a g(T — t)-isometry we have: 

((df(T - t, d t {g{T ~ t)) #T -We 4 ))) l=1 .. n , (t^^Wj^R- 

= (Zt 9t(g(T - mufe^ v r -V(T - 1, .)K (^)~ 1 w^) R . 

= <Ei - t)Wle t , V T -*f{T - t, .))Ufe h W^) T _ t 

= (d t (g(T ~ t))# T ~\Wl t v), V^fiT - t, .)> r -t, 
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Consequently: 

d(df(T-t,.) X T {x) (Wl t v)) 
= -j t df(T-t : .)(Wl t v)dt 

- \ <V T -7(T -*,.), d t (g(T - t))#^(W^)) T _ 4 dt 
+ ^A T _ t df(T-t,.)(Wl t v) dt 

- l -({df{T - t, U?)) t {Uj)~\mc g(T _ t) -d t (g(T - t)))#^-')(Wj>dt)j 



-d/(T - t, -)(W^) dt + -A T _ t df(T - t, .)(W^) dt 



1 



-^/(T-^Ricf^-^W^^dt. 



But recall that / is a solution of: 



W ~ 2 A</ ' 



so that 

-^-df(T-t,.) = ~dA T _ t f(T-t,.). 

We shall use the Hodge-de Rham Laplacian Qr_t = —{d5x~t + Sr-td) which com- 
mutes with the de Rham differential, and we shall use the well-known Weitzenbock 
formula (pH EJ), which says that for 9 a 1-form, U T _ t 9 = A T _ t - Ric T _ t 9. We 
get: 

dA T _ t f(T — t, .) = dD T _ t /(T-t,.) 

= □ T _ t d/(T-t,.) 

= A T _ t d/(T-t,.)-RicT- t d/(T-t,.). 

Finally: 

d(df(T-t,.)xf(,)(Wj^)) "= |Ric T _ t d/(T-t,.)(W^)dt 



f o, 2 



1 (v t -V(t - 1, .), mcjFr*(w3>)>T-t ^ 



by duality; for a 1-form and for v G TM: 

Ric(0)(u) = Ric(fl # ,^ 

where {9*,v) = 9{y) . □ 
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Remark : For the forward Ricci flow, we have: 

//lt*d((//l t y 1 wl) = o. 

For the backward Ricci flow, we have: 

lilt * ddl/lr^lt) = - Ric#^(Wj t ) dt. 

When the family of metrics is constant, we have the usual damped parallel 
transport, wich satifies: 

l/o,t * d((//o,0 _1 Wo,t) = ~\ Ric # (W , t ) dt. 



Remark : Roughly speaking, the result says that the deformation of the met- 
ric under Ricci flow makes the damped parallel transport behaves like the damped 
parallel transport in the case of a constant metric with flat Ricci curvature. 

For the heat equation under the forward Ricci flow, we take the probabilistic 
convention: 

r d t f(t,x) = ±A t f(t,x) 

{ &«.,■ = -Ricy (3.3) 
( f(0,x) = f (x) 

We shall give a Bismut type formula and a gradient estimate formula for the 
above equation. For notation, let T c be the maximal life time of the forward Ricci 
flow g(t)te[o,T c [, solution of (12.21) . For T < T c , Xj is a g{T — t)-Brownian motion 
and //qj the associated parallel transport. In this case, for a solution f(t, .) of 
(13.31) . f(T — t,Xf(x)) is a local martingale for any x G M. When going back in 
time, one has to remember all deformations of the geometry. 

We now recall a well known lemma giving a Bismut type formula (e.g. [8]). 
Let f(t, .) and g{t) be solution of (13.31) . T < T c , and Xj{x) a g(T — t)-Brownian 
motion. 

Lemma 3.3 For all W 1 -valued process k such that k e Lf oc (W) where W is some 
M. n -valued Brownian motion, and for all v G T X M , 

N t = df{T-t,.) X T {x) {Uf)[{U^r l v-jlk T dr} 
+ f(T-t,X?(x))j*(k r ,dW) Rn 

is a local martingale. 
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Proof: The first remark after theorem 13.21 yield that the first term is a semi- 
martingale. By ltd calculus we get: 

d(f(T - t, Xj(x))) = df{T - t, ) xnx) U t e l dW i . 
With (Zj)j=i„ n a g(T)-orthonormal frame of T X M, we write N t as: 

Nt = OT(r-t 1 ^ w (TO- 1 )«(^-J?(TO,yT*) 

+ f(T-t,Xf(x))^(kr,dW)^ 

with[32J 

dN t "= E i (df(T-t 7 .) xnx) (U^(U^r%)(-(U^k t )J l ) T dt) 
+d{f(T-t,X?(x)))(kt,dW)vn 

+ df(T - 1, .) xT{x) {Ufl l )dW\Y, 3 KdW?) 
0. 

□ 

Remark: Since T is smaller than the explosion time T c , and by the com- 
pactness of M, N t is clearly a true martingale, so we could use the martingale 
property for global estimate, or the Doob optional sampling theorem for local es- 
timate (e.g. [23]). 



Corollary 3.4 Let v G T X M , and take for example k r = — -l[o,T](r) then: 

d/(r ) .)^ = ij]E[/ (^(x))((^)- 1 ^,e i ) H «Wi(r)]. 

i 

Proof: With the above remark, N t is a martingale. The choice of k r gives 

(Uq)~ x v — J Q T k r dr = 0; the result follows by taking expectation at time and T. 
□ 

We can give the following estimate for the gradient of the solution of ( 13.3( 1: 
Corollary 3.5 Let ||/||oo = sup M |/ |. For T < T c : 

sup || V T /(T, x)\\t is decreasing in time 
xeM 

and: 



sup ||V T /(T,*)|| T < 

x£M V J- 
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Proof : Take x G M such that || V T /(T, x) \\t is maximal. Using the damped 
parallel transport 13.21 we obtain that for all v G T^M: 

d/(T-t,X f T (x))W^, 

is a local martingale. By compactness, this is a true martingale. Taking v = 
V T /(T, x) and averaging the previous martingale at time and t we get: 

|| V T /(T,x) |||=E[(V T -V(T-^X t T (x)),Wj^) T - t ]. 



Using l3~2l we get the first result. 



If we choose k r 



l[o,T]( r ) in 13.31 then N t is a martingale. Taking 



expectations at times and T, we obtain 



1 



df(T,.) x v = -E[f (X±(x)) I (Ut)-\dW) Rn }. 



For x G M and = V /(T, x), Schwartz inequality gives 



V7(T,i) ||^< 



T 

(f/ o r )^ 1 t;,^) i 



We have: 



E 



T 

{UT)- l v,dW)^ 



T || v 112 



T • 



The result follows. □ 



For geometric interpretation, let us give an example of normalized Ricci flow 
for surfaces (which is completely understood e.g. [5]). We are interested in this 
example because the equation for the scalar curvature under this flow is a reaction- 
diffusion equation which is quite similar to the heat equation under Ricci flow. We 
will give a gradient estimate formula for the scalar curvature under normalized 
Ricci flow which gives in the case x(M) < (the easiest case) the convergence of 
the metric to a metric of constant curvature. 

The normalized Ricci flow of surfaces comes from normalizing the metric by 
some time dependent function to preserve the volume. Let M be a 2-dimensional 
manifold, R(t) the scalar curvature, r = J M Rtd/i t / ')Ut(M) its average (which will 
be constant in time, as topological constant, e.g. Gauss-Bonnet). We get the 
following equation for normalized Ricci flow: 
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Remark : Hamilton gives a proof of the existence of solutions to this equa- 
tion, defined for all time ( e.g. [5])). 

Recall that (e.g. [5]) the equation for the scalar curvature R is: 

^R = A t R + R(R- r). 

Proposition 3.6 Let T e M, Xf(x) be a \g{T — t)-BM(x) , //^ t the parallel trans- 
port,, v G T X M and ip t v the solution of the following equation: 

//ltd ((//j t rw) = -(l r - 2R ( T - ^ T (*))) ftvdt 

ipo = Id^M • 

Then dR(T — t, ■)x T {x) l Pt v is a martingale and: 

\\V T R(T,x)\\ T < sup ||V° J R(0,x)|| e-f rT E[e^ T2iJ ( T -^ T W) rf *]. (3.4) 

M 

Proof: The proof is similar to the one in 13.21 the difference is the reaction term: 
R(R—r). For notations and some details see the proof of l3T2l Take F:xh-> x(x— r), 
then: 

^R = A t R + F(R). 

We write: 

dR(T - t, .) \ x t [x) <p t v = (d~R(T - t, Uj), {Ujy l Vt v) W 2 
where Uf is a diffusion on T{M) with generator 

A%_ t + ±( r -R(T- t, n.))g(T - t)(ev ei ., ev ej .)V hJ {.). 

Using theorem 13.21 we have: 

d{dR(T -t,C^),(C^)-Vtv)R a 
= (d(dR(T - t, Ul)\ {U?)- l <ptv)w 
+ (dR(T - t, Uj), d{(V*)- x tp t v))w 

d = \—(dR(T - t, .)) + A T . t dR{t -t,.) + -(r- R(T - t, n.))dR(T - t, .)] (<p t v) dt 
+ (dR(T -t,Uj),d{(Ujy l ip t v)) R , 

Using Weitzenbock formula and the equation for R we get: 

^-dR(T-t, .) = -[A T - t dR(T-t, .)-Ric T _ t dR(T-t, .)+F'(R(T-t, .))dR(T-t, .)] 
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Recall that for the surface: 

Ric T _ t dR(T - t, .) = -R(T -t,.)dR(T -t,.), 

consequently 

d(dR(T -t,[/n,(CrW 

^ (I r - F'(i?(T - t, .))dR{T - t, .))(**;) eft + (d~R(T - t, Uj),d{{Uj)- l ^ t v)) R , 
d = (ir - 2i?(T -*,.) + r)di?(T - t, .))(M d* 
+ (dR(T - i, Uj), (U?)-\-~r + 2i?(T - t, 

where we used the equation of <^t> in the last step. 

For the second part of the proposition, with the equation for <p t v we have: 

(p t v \\ 2 T _ t ) = {4R(T -t,Xj(x) -3r) || ip t v \\ 2 T _ t dt, 



so that 



VTV \\l=\\ if Q V \\ 2 T e -3rT e jTm(T-s,XT( X ))ds_ 



Take v = VtR{T, x) and average at time and T (it is a true martingale 
because all coefficients are bounded) to get: 

\\V T R{T,x)\\ T < sup ||V J R(0,x)||oe"f rT E[e/o T2 ^-^ T (-))^]. 

M 

□ 

Remark: For reaction-diffusion equations we can find by this calculation 
the correction to the parallel transport leading to a Bismut type formula for the 
gradient of the equation: 

lf = A t f + F(f), (3.5) 

where A t is a Laplace Beltrami operator associated to a family of metrics g(t). 
Let Xf(x) be a \g(T — t) — BM(x), //^ t the associated parallel transport and 
v G T X M. Consider the covariant equation: 

//lA/flr^tV = - ( Ric # < T -' -\ [^(g(T - t))] *' T ~ l - F\f)) Q t v dt 
Then for / a solution of f!3.5p and v 6 T X M we obtain that: 

df(T-t,.)e t v 

is a local martingale. 
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Corollary 3.7 For x(M) < 0, there exists C > depending only on g(0), such 
that: 

|| W T R(T, x)\\ T < sup \\V°R(0, x) || e^e 2 ^ 2 ^. 

Proof: We use proposition 5.18 in [5]. In this case we have r < and a con- 
stant C > depending only on the initial metric such that R(t, .) < r + Ce rt and 
the estimate follows from previous proposition. □ 

Remark: For the case x{M) < we obtained an estimate which decreases 
exponentially. For the case x(M) > one could control the expectation in (13.41) . 



4 The point of view of the stochastic flow 

Let g(t)[o,T c [ be a C 1 ' 2 family of metrics, and consider the heat equation: 

\ /(0,x) = / (x), 

where fo is a function over M. We suppose that the solution of this equation exists 
until T c . For T < T c , let Xf be a g(T — t)-Brownian motion and //g t the associated 
parallel transport. 

We will build (c.f. (14. 21) ) a family of semimartingales (T — t,X~[(x)) such 
as Xj{x) is a g(T — t)-BM(i) for all x nearby x and such that the family of 
martingales f(T — t,Xf{x)) x is differentiable at xq with respect to the parameter 
x. However, in this section, we will not do it directly using stochastic flows in the 
sense of [20]. Instead, we will use differentiation of families of martingales defined 
as limit in some semi-martingale space (the topology is as in [TT] which has been 
extended by Arnaudon, Thalmaier to the manifold case [1], [3], [I], [2]). 

We work in the space-time IxM, its tangent bundle being identified to TIxTM 
endowed with the cross connection V = V ® Vr-t where V is the flat connection. 
Let Xf(xo) be a g{T — t)-BM started at Xo, and define Y t (xo) = (t,X"[(xo)) a 
/ x M- valued semi martingale. From now on P][ Y stands for the parallel transport 
along the shortest V-geodesic between nearby points X G / x M and Y G / x M 
for the connection V. 

Let c a curve in I x M, we write P s v for the V parallel transport along c and 
for a curve c in M we denote by the V T_S parallel transport along c. We 

also denote it : I x M — > M the natural projection. 

For a curve 7 : t — > (s,x t ) in I x M, where s is a fixed time, we have the 
following observation: 
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Define the Ito stochastic equation in the sense of [1 3j : 

cFY t (x) = P? t{xo)Mx) d*Y t (x ) (4.2) 

Remark: The above equation is well defined, for x sufficiently close to xo, 
because dT-t(X t (x), X t (x )) is a finite variation process, with bounded derivative 
(by a short computation and [TS], [6]). 

Let // ot be the parallel transport, associated to the connection V, over the semi 
martingale Y t (xo). 

In the next lemma, we will explain the relationship between the two parallel 
transport // ot and //^ t . 

Lemma 4.1 Let (ei)i=i.. n be a orthonormale of (T X0 M, g(T)) then 

diU/lr'd-K//^ ei) = \{lll t r\j9{T - t))* T ~\dnj/^ *)) dt. 
Proof: The parallel transport // 0i does not modify the time vector, i.e., 

//rio(°»-) = (o.-), 

as can be shown for every curves, and hence for the semi-martingale Y t by the 
transfer principle. 

We identify T = {(0, v) G T {0jXo) I x M} and T X0 M with the help of (0, v) i — ► v. 
Hence 

U/ltY^JUt : f -> T X0 M 

becomes an element in A^ n , n (M). 

Recall that //£ t = UjU^'' 1 . By definition of D 5 '* given in (j 1 . 2f) . We get using 
the shorthand = L/Jej, with (ej)j = i.. n an orthonormal frame of R n , 

*d(U/Z t )- 1 dn// J = *d(((//l t )- 1 dn// 0tt e i ,e j ) T ) i , j 
= *d((dTT// je i J/l t e j ) T „ t ) iJ 
= ((D s ' T - t dn// ^U?e J ) T _ t 

+ ^(g(T-t))(d 7 r// 0>t e i ,U?e j )dt 

+(dn// Q!t e i ,D s ' T - t U^e j ) T _ t ) . 
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We also have: 



Where we have used in the last equality the fact that // 1 is the V horizontal lift 

of Y t . The third one may be seen for curves, it comes from the definition of V. 
Following computations similar to one in the first section, we have by (11.51) : 

*d(u/o,t)~ 1 dW/ 0>t )ij = yr-tirt^.^)* 

= ^{T-t^ditJ/^Ufe^dt 

+ (dir//^, -\ Eti W t 9( T - Ul^Uje^T-t dt 

= ^g(T-t)(dir// 0it e i ,UTe j )dt 

~ \ ELi b{T - t)(U?e jt U?e a )(d*// 0tt e it Uj~e a ) T . t dt 

= l^giT-t^J/^Ufe^dt. 

In the general case, and by previous identification: 

d((//l t r 1 d7r//^ t )(0,e i ) = lj2^9(T-t)(d7r// 0tt e i ,U^e j )e j dt (4.3) 

3 

= \{/llt)-\j9{T - t))* T -\dv//^ e,)) dt. (4.4) 

□ 

Differentiating (14.21) along a geodesic curve beginning at (0, xo) with velocity 
(a, v) G T I x T X0 M and using corollary 3.17 in [3] we get: 

IUAlUl TY tM) = ~R{TY t {a,v),dY t (x ))dY t (x ), 

where R is the curvature tensor. 
Let v G T X M we write: 

TX t v := dnTY t (0,v). 
In a more canonical way than theorem 13.21 we have the following proposition. 



V 

( 

0. 



(ddn(dev e A*d// ,)Y) 
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Proposition 4.2 For all v £ T X M we have: 

diU/lr'TXtv) = \{//l t T\{j9{T - t)) - mc T _ t )#> T -\TX t v) dt. 

Proof: For a triple of tangent vectors (L t , L), (A t , A), (Z t , Z) £ TI x TM, we 
have: 

R((L t , L), (A u A))(Z t , Z) = (0, R T -t(L, A)Z). 

Hence, according to the relation dY(xo) = (dt, *dX t ) = (dt, ll% t e% * dW l ) and the 
definition of the Ricci tensor: 

//oA//o!tTYt(0,v)) = -^(0,Ric #T -'(TX^))rft. (4.5) 

In order to compute in R n , we write: 

U/lr'TXtV = ((//^)- 1 ^// 0ii )(// >r i (0, «)). (4.6) 

By (|4.5p . we have rf(// t TY t (0, t>)) £ dA where A is the space of finite variation 
processes. We get: 

d((//lt)- 1 TX t v)=d((//lA^dn//^ 
By (14.61) and lemma ED we get: 

d({//l t )-'TX t v) = *d{U^- 1 dn//^{//^ t TY t (0,v)) 

+ (u/hr^j/^^du/^TY^v)) 

= ^((//^J" 1 ^// 0i j(//ojry,(o^)) 

- |((//ot) _1 ^)(0, Ric # '^*(TX^) dt 

- W/l^Kic^iTX^dt. 

□ 

For all /o £ C°°(M), and for f(t,.) a solution of dHHJ), /(T - t,Xj(x)) is a 
martingale, where (T — t, Xj(x)) = Y t (x) is built as in (I4.2p . We have the following 
corollary which agrees with theorem 13.21 

Corollary 4.3 For all v £ T X M: 

df(T - t, Xj(x))v = df(T - t, .)xT ( »)//o>» 

is a martingale. 
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Proof: By differentiation under x of f(T — t, X~[(x)), we get a local martingale. 
According to [3] and by chain rule for differential we get the corollary. This result 
matches 13.21 after using the above proposition. □ 

In an canonical way, we have the following result. 
Theorem 4.4 The following conditions are equivalent for a family g{t) of metrics: 
i) g(t) evolves under the forward Ricci flow, 
ii) For all T < T c we have //^ t = W^ t = TX t . 

Hi) For all T < T c , the damped parallel transport W^ t is an isometry. 

Proof: Here, the forward Ricci flow has probabilistic convention (12 .21) . The 
result follows by the equation of g{t) and by proposition 14.21 and theorem 13.21 □ 



5 Second derivative of the stochastic flow 

We take the differential of the stochastic flow in order to obtain a intrinsic mar- 
tingale. We take the same notation as the previous section, and g{t) is a family of 
metrics coming from a forward Ricci flow. Let Xf(x) be the g(T — t)-BM started at 
x, constructed as in the previous section by the parallel coupling of a g(T — t)-BM 
started at x ( 14.21 ). V and Y t (x) = (t,Xf(x)) as before, define the intrinsic trace 
(that do not depend on the choice of Ei as below): 

TrV.TX t (x )(.) :=rf7r^V ( o, e ^(x)(0,^(x))-TFi(x)V ( o, ei )(0,^(x))') 

where (e^) is a (T Xo M, g{T)) orthonormal basis, Ei are vectors fields in YTM such 
that Ei(x ) = ei and V (o )ei )TY t (x)(0, Ei(x)) is a derivative of a bundle-valued 
semi-martingale in the sense of ([1], [3], [1]). By 14.41 : 

Tr VTX t (x )(.) := dn £ V^)TY t {x)(0, E,{x)) - V (0 , ei) (0, E^x))) 

i i 

Theorem 5.1 Let L t := Tr VTX t {x ){.) be a (T Xo M, g(T)) -valued pro- 

cess, started at 0. Then: 

i) L t is a (T X0 M, g(T))-valued martingale, independent of the choice of Ei. 

ii) The g(T)- quadratic variation of L is given by d[L,L] t =|| Ric T ~*(X t (a;o)) Wr-t 
dt. 
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Proof: Recall that by the same construction of the previous section: 

D(TY t (x)(p,Ei(x))) = -±R(TY t (x)(O t Ei(x)),dY t (x))dY t (x). 

By the general commutation formula (e.g. theorem 4.5 in [4]), and by the previous 
equation which cancels two terms in this formula, we get: 

DV {0 , ei) (TY t (x)(0,E t (x))) =V ( o, ei)J D(TF t (x)(0,^(x))) 

+ R(o?Y t {x ),TY t (x )(0,e i ))TY t (x Q ){0,e i ) 

- ^R(dY t (x o ),TY t (x )(0, e*), dY t (x o ))TY t (x )(0, e*) 

= - ^v ( o,e 4 )(^(ry t (x)(o,^(x)),dy t (x))dy t (x)) 

+ R(d*Y t (xo),TY t (x Q )(0,e i ))m(xo)(0,e i ) 

Taking trace in the previous equation we can go one step further. Recall that 
(ei)i=i..n is a orthogonal basis of (T XQ M, g(T)), and write for notation: 

mcl x) (v) = (o,mc* T - t (dnV)), 

then: 

E< 5v M (3y t (i)(o^(i))) 

= ~\ Ei V^e^RicJ^^F^x)^^))) 

+ ^(rf v y(x ), TF (x )(0, e,))TF 4 (x )(0, e<) 
-|E i (Vdy i( , )^)(TF t (x )(0,e J ),dy(xo))TF t (x )(0,e i ) 

= -| Ei(V(Ty t (xo)(o, ei ))Ric # )(Tr t (xo)(0, e*)) 

-KR i 4(,o)(E,V { o,e l )TF 4 (x)(0,^(x)))) + Ric5 (xo) (^(xo)) 
-lEi^^i?)^^)^,^),^^))^^)^,^). 

In the last equality, we use the chain derivative formula, and derivation is taking 
with respect to x . We will make an independent computation for the last term in 
the previous equation. Let Tr stand for the usual trace: 

E t (^dY t ( xo) R)(TY t (x )(0,e l ) J dY t (x ))TY t (x )(0,e l ) 

= Ei(0, (V^ J R T -')(TX 4 (x )e,,ciX t )TX t (x )e i ) 

= E M (0, (V^ e .i? T -*)(rX t (xo) ei , //l t e,)TX t (x ) ei ) dt 

= E,(0, Trx^V^^^C//^-)) dt 

= E,(0, (V*^ Tr^i^) (//£>,■)) dt 

= - E,(o, (vg ej Ric# T -*)(//2>i)) di, 
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where we have used in the second equality the fact that in case of the forward 
Ricci flow //q t is a g(T — t) isometry and dX = //^ t ejdWK In the last equality we 
use the commutation between trace and covariant derivative (for example [2T], or 
p[9]). Note that: 

E l (V(Ty t (xo)(o, ei ))R ic# )( Tr t( a; o)(0,e i )) 
= E*(0, iy T TX%)e, Ricf^TX^e,)) dt 

Hence, using WM 

(o,e t )TY t (x)(0, Ei(x))) 
= -|(Ri(Vt(*o)(Z)i ^(o,ei) TY t(x)(Q, E i(x)))) + Ric* t{xo) (d7Y t (x )) 

Write, for simplicity, B for ^V^^^^^^)). We compute: 

di/ft-'dnB) = d([//l :t ' 1 d7r// 0tt ][(//o,r 1 B]) 

= \//l;\d t g(T -t))*- T -\d*B)dt 

+lf$\-\d*{m?{B)) + dn(Ric* {xo) (d*Y t (x )))) 
= //l-\dnRic*(d^Y t (x ))) 
= E i //Jr 1 Ric$J ) (//^e i )^, 

where we have used lemma 14.11 in the first equality. We get a intrinsic martingale 
that does not depend on Ei, starting at 0. By the definition in theorem 15.11 and by 
the formula preceding theorem 15.11 the above calculations yield: 

L * = f E Hit 1 Ricf^U/l^dW* - d7r(J2 V ( o, e! )(0, Ek(x))). 

For the <?(T)-quadratic variation of L t we use the isometry property of the parallel 
transport; we compute the quadratic variation: 

d[L,L] t = (//^Ric*^ 

= Ei II Ri C xT(x)dlo,t e i) lls(T-t) dt 

- Ill "Ric #T_ * III 2 dt- 

- Ill ruL xT(x) WlT-t at ' 

where |||.||| is the usual Hilbert-Schmidt norm of linear operator. By the indepen- 
dence of the choice of the orthonormal basis we can express this norm in terms of 
the eigenvalues of the Ricci operator: 

d[L,L} t = J2^(T-t,X?(x))dt. 
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□ 

Remark: We could choose Ei such that V(o, ei )(0, Ei(x)) = that do not 
change the martingale L, but give a simple version. 

Remark : This martingale can be used to look at the behavior at point where 
the first singularity of the Ricci flow occurs, i.e. where the norm of the Riemannian 
curvature explodes ([5], [2]). 
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